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Abstract. We examine the spatial distribution of brown dwarfs produced by the decay of small- N stellar systems as expected 
from the embryo ejection scenario. We model a cluster of several hundred stars grouped into 'cores' of a few stars/brown 
dwarfs. These cores decay, preferentially ejecting their lowest-mass members. Brown dwarfs are found to have a wider spatial 
distribution than stars, however once the effects of limited survey areas and unresolved binaries are taken into account it can 
be difficult to distinguish between clusters with many or no ejections. A large difference between the distributions probably 
indicates that ejections have occurred, however similar distributions sometimes arise even with ejections. Thus the spatial 
distribution of brown dwarfs is not necessarily a good discriminator between ejection and non-ejection scenarios. 
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1. Introduction 



Brown dwarfs are observed to constitute some 15-25% of 
the objects in young star forming regions (e.g. Bricefio et 
al. 2002; Muench et al. 2003; Luhman et al. 2003; Luhman 
2005). However, the formation mechanism(s) of brown dwarfs 
are currently unclear. The two most popular models are that 
(a) brown dwarfs form like stars but in very low-mass cores 
(e.g. Greaves this volume; Padoan this volume), or (b) brown 
dwarfs are stellar embryos that are ejected before than can 
accrete sufficient mass to become stars (Reipurth & Clarke 
2001,2003). (Also see Whitworth & Goodwin (2005) and 
Whitworth & Goodwin (this volume) for a review of other 
possible mechanisms). 

Probably the most popular model for brown dwarf forma- 
tion (at least among theorists) is the ejection scenario (Reipurth 
& Clarke 2001). This model appears to be supported by sim- 
ulations of star formation in turbulent cores (e.g. Bate et al. 
2002, 2003; Delgado Donate et al. 2004; Goodwin et al. 2004a,b)2. Method 

The ejection scenario has three significant problems. Firstly, 



recently become more acute with evidence that the brown 
dwarf binary fraction may be significantly higher than pre- 
viously thought (Pinfield et al. 2003; Jeffries & Maxted this 
volume). Finally, brown dwarfs are ejected with velocities of 
order 1 km s _1 and it may be expected that brown dwarfs are 
more widely distributed than stars, in contrast to observations 
(Bricefio et al. 2002; see also Luhman 2005). 

In this paper we concentrate on the final problem - the 
spatial distribution of brown dwarfs. In the ejection scenario, 
cores form several objects in an unstable high-order multi- 
ple system. Dynamical interactions typically eject the lowest- 
mass members of the system until a stable binary or hierar- 
chical multiple remains. We perform A^-body simulations of 
a cluster of decaying small-iV 'cores' to examine if the distri- 
butions of brown dwarfs and stars are different or if the veloc- 
ity dispersion between cores effectively conceals the ejection 
of brown dwarfs and low-mass stars. 



it is not clear if ejected brown dwarfs are able to retain the 
significant discs which are observed around at least some 
young brown dwarfs (e.g. Jayawardhana et al. 2003). Sec- 
ondly, whilst ejections do produce some brown dwarf -brown 
dwarf binary systems (e.g. Bate et al. 2002) a large popu- 
lation of brown dwarf -brown dwarf binaries would be dif- 
ficult to explain via the ejection scenario (this problem has 



We simulate the A^-body evolution of a cluster containing 
many 'cores' containing several stars using the NBODY6 code 
(Aarseth 2003). 

A core is a small-A*" system of stars and brown dwarfs 
in which we assume accretion has finished and the gas reser- 
voirs have been exhausted before dynamical evolution occurs. 
This is a simplification to avoid dealing with the complex gas 



3 RESULTS 



dynamics of the combined accretion/ejection phase (e.g. Del- 
gado Donate et al. 2004; Goodwin et al. 2004a,b; see Umbreit 
et al. 2005 for an TV-body treatment of this problem). How- 
ever, we believe that it retains the key physics of the current 
problem - the ejection of low-mass members of cores which 
are moving relative to one-another within a cluster. 

A cluster is initially composed of A coro cores each con- 
taining stars and/or brown dwarfs (so that the total num- 
ber of stars in the cluster is TVtot = A cor c x A»). For a typical 
distribution the inter-core velocity dispersion is around 1 km 
s" 1 (of order the ejection velocities of low-mass stars and 
brown dwarfs - see below). 

Within a core the masses of stars (for brevity 'stars' gen- 
erally means both stars and brown dwarfs) are randomly sam- 
pled from a Kroupa (2002) IMF of the form 



N(M) 




0.02 < M/Mq < 0.08 
0.08 < M/Mq < 0.5 
0.5 < M/Mq < 5 



The N* stars are distributed randomly within a region of ra- 
dius i? CO rc and given random velocities which are scaled such 
that the core is in virial equilibrium. We note that a random 
sampling of the IMF in this way may not reproduce the cor- 
rect (post-ejection) multiple system properties (Kroupa this 
volume). 

We choose R COIC = 200 au as the typical scale on which 
stars are expected to form (e.g. Goodwin & Kroupa 2005). 
This is the scale at which a collapsing core will reach the 
critical density at which the minimum mass for fragmenta- 
tion is reached (~ 10~ 13 g cm" 3 ). We note that such a length 
scale has an observational basis as the peak of the T Tauri bi- 
nary separation distribution occurs at ~ 100 au (e.g. Mathieu 
1994; Patience et al. 2002). 

The cores so established are expected to decay on a timescale 
of < 0.1 Myr rapidly if A sta r > 2 by ejecting the lowest- 
mass members of the system until a stable hierarchical sys- 
tem or a binary is formed (Anosova 1986; see also Sterzik 
& Durisen 2003; Hubber & Whitworth 2005; Goodwin & 
Kroupa 2005; Umbreit et al. 2005). We ignore the interaction 
of the protostars and the ambient gas, as all we are interested 
in is the post-ejection velocity distribution and we find that 
the ejection velocities of low-mass stars and brown dwarfs 
are of order 1 km s _1 , similar to that found in more detailed 
simulations including gas (e.g. Umbreit et al. 2005). 

Cores are then placed within the star cluster by position- 
ing them within a virialised Plummer sphere with a virial 
radius of 1 pc following the prescription of Aarseth et al. 
(1974). 



3. Results 

We analyse the relative distributions of brown dwarfs and 
stars by comparing the distances to the nearest neighbours 
(following Brinceno et al. 2003); using the mean distances to 
the nearest neighbour (nearest neighbour distance, or NND). 

There are two important observational biases that must 
be included when analysing the data. The first is unresolved 



binaries, which we include by ignoring the secondary com- 
ponent of any system if it is within 250 au of the primary 
(roughly 2 arcsec at the distance of Taurus). The second is 
that surveys are over limited areas: cluster members at great 
distances from the cluster centre will normally not be found 
in surveys, and - even if they are - would be difficult to un- 
ambiguously relate to the cluster without additional proper 
motion studies. Therefore we restrict ourselves to stars and 
brown dwarfs within a projected distance of 5 pc from the 
cluster centre (this corresponds to an area of 13 square de- 
grees at the distance of Taurus). 

In Fig. 1 we compare the stellar and brown dwarf NNDs 
for clusters with a total of N to t — 400 stars and = 1 
(ie. no decay of groups within cores) and A* = 4 (generally 
decay into a binary and two single stars). 

When A* — 1, the NNDs of the stars and brown dwarfs 
are indistinguishable, over the 20 Myr of cluster evolution 
that is followed, the NNDs increase somewhat as the cluster 
expands slightly through 2-body interactions. Over a signif- 
icant time, we would expect the brown dwarf NND to become 
larger than the stellar NND as low-mass stars and brown dwarfs 
gain a higher velocity dispersion as equipartition is estab- 
lished through 2-body encounters (however, in reality the clus- 
ter will probably disperse long before this becomes impor- 
tant). 

In contrast, when A* = 4, the effect of the very rapid 
decay of the small- A cores is to produce a population of 
more widely dispersed brown dwarfs and low-mass stars with 
a larger NND. (It should be noted that not accounting for un- 
resolved binaries and limited survey areas makes these dif- 
ferences significantly more extreme). Towards the end of the 
simulation the NND for brown dwarfs is seen to drop signifi- 
cantly. This is caused by a number of brown dwarfs escaping 
from the cluster at late times, either because they were ejected 
with an initially low velocity or have gained velocity due to 
later encounters in the cluster. 

Thus it would appear that the ejection of brown dwarfs 
from small- A cores can produce a significant and observable 
effect in the spatial distributions. However, ejections do not 
always produce a significant difference in the spatial distri- 
butions. In Fig. 2 we show four more simulations with A* = 
4 (making 5 in total including the simulation from Fig. 1). 
The only difference between these simulations is the random 
number seed used to generate the initial conditions. In the fi- 
nal simulation in particular there is no significant difference 
between brown dwarfs and stars at any time. 

The most extreme differences between the NNDs of brown 
dwarfs and stars occurs at 5 - 10 Myr: after ejections have had 
time to significantly disperse the brown dwarfs, but before 
most of them have escaped the inner regions of the cluster 
(travelling at 1 km s _1 this should take ~ 5 Myr). However, 
most young clusters are observed at ages of ~ 1 Myr (e.g. the 
Taurus observations of Briceno et al. 2002) at which point 
only 1 of the 5 simulations shows a very significant differ- 
ence. 
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4. Conclusions Umbreit, S., Burkert, A., Henning, T., Mikkola, S. & Spurzem, R. 

2005, ApJ, in press lastro-ph/0501075 l 

When ejections are not important (iV* = 1) the spatial distri- 
butions of stars and brown dwarfs are (unsurprisingly) very 
similar. When ejections become important (iV* = 4) then the 
spatial distributions can show significant differences. How- 
ever, these differences can disappear altogether in some clus- 
ters depending on the exact details of the initial clustering. In 
addition, at the young age of most well-studied young clus- 
ters a significant difference in the spacial distributions is seen 
in only 1 out of 5 simulations. 

Thus if brown dwarfs and stars have different spatial dis- 
tributions it is probably a signature of ejections, however the 
lack of a difference does not necessarily exclude the ejection 
scenario. 

We suggest that binarity is a far stronger discriminator 
between models since (as yet) there is no way in which to 
make significant numbers of brown dwarf -brown dwarf bina- 
ries from ejections. However, we are currently investigating 
if a more sophisticated statistical analysis of positions may 
yield a more robust discriminator. 

The discovery of many more brown dwarfs in Taurus ex- 
tending over a wide area (Guieu et al. this volume) does sug- 
gest that ejections may have been responsible for at least 
some of the Taurus brown dwarfs. 
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Fig. 1. The average nearest neighbour distances (NND) of stars (dahsed-line) and brown dwarfs (full line) for clusters with 
A/tot = 400 with JV* = 1 (left panel) and N* = 4 (right panel). 
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Fig. 2. The average nearest neighbour distances (NND) of stars (dashed-line) and brown dwarfs (full line) for clusters with 
Ntot — 400 and iV* = 4 for four different realisations (different random number seeds). 



